Cosmological Observations: Averaging on the Null 

Cone 



A. A. Coley 

Department of Mathematics and Statistics 
Dalhousie University, Halifax, NS B3H 3J5, Canada 
aac@mathstat. dal. ca 



Abstract 



The universe is not isotropic or spatially homogeneous on local scales. The aver- 
aging of local inhomogeneities in general relativity can lead to significant dynamical 
effects on the evolution of the universe and on the interpretation of cosmological data. 
In particular, all deductions about cosmology are based on light paths; averaging can 
have an important effect on photon propagation and hence cosmological observations. 
It would be desirable to describe the physical effects of averaging in terms of observa- 
tional quantities and focussing on the behaviour of light. Data (e.g., matter terms, such 
as the density or galaxy number counts, which are already expressed as averaged quan- 
tities) is given on the null cone. Therefore, it is observationally meaningful to consider 
light-cone averages of quantities. In principle, we wish to describe the cosmological 
equations on the null cone, and hence we need to construct the averaged geometry on 
the null cone. However, we argue that it is still necessary to average the full Einstein 
field equations to obtain suitably averaged equations on the null cone. Since it is not 
the geometry per se that appears in the observational relations, we discuss whether it 
is possible to covariantly 'average' just a subset of the evolution equations on the null 
cone, focussing on relevant observational quantities. We present an averaged version of 
the scalar null Raychaudhuri equation, which may be a useful first step in this regard. 
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1 Introduction 



Cosmological observations [H [2], based on the assumption of a spatially homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) model plus small perturbations, 
are usually interpreted as implying that there exists dark energy, the spatial geometry is 
flat, and that there is currently an accelerated expansion giving rise to the so-called ACDM- 
concordance model. Although the concordance model is quite remarkable (at least if the 
idea of dark energy can be tolerated), it does not convincingly fit all data. Essentially, 
structure is more evolved on large scales than predicted by the ACDM-concordance model. 
There is strong evidence for coherent bulk flows on scales out to at least 300/i _1 Mpc [3] 
and there is observational evidence of anisotropy in the Hubble expansion rate, suggesting 
large scale peculiar motions [I]. In addition, there are correlations between galaxy surveys 
and the cosmic microwave background (CMB) data [5], overestimates of the amplitude of 
the matter power spectrum [6] and non-Gaussianities in the CMB [TJ. Unfortunately, if the 
underlying cosmological model is not a perturbation of an exact flat FLRW solution, the 
conventional data analysis and their interpretation is not necessarily valid. For example, the 
standard analysis of type la supernovae (SNIa) and CMB data in FLRW models cannot be 
applied directly when backreaction effects are present, because of the different behaviour of 
the spatial curvature [5]. 

Supernovae data dynamically requires an accelerating universe. However, this only im- 
plies the existence of dark energy if the universe is (approximately) FLRW. Thus the isotropic 
and spatially homogeneous ACDM model is a good phenomenological fit to the real inho- 
mogeneous universe, as far as observational determinations of the expansion history of the 
universe [9]. However, this does not imply that a primary source of dark energy exists, 
but only that it exists as far as the phenomenological fit is concerned. Supernovae data 
can be explained without dark energy in inhomogeneous models, where the full effects of 
general relativity (GR) come into play. The apparent acceleration of the universe is thus 
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not caused by repulsive gravity due to dark energy, but rather is a dynamical result of inho- 
mogeneities, either in an exact solution or via averaging effects (due to the back-reaction of 
inhomogeneities). Indeed, it has been indeed shown that the Lemaitre-Tolman-Bondi (LTB) 
solution can be used to fit the observed data without the need of dark energy [10] , although 
it is necessary to place the observer at the center of a rather large-scale underdensity. 

Therefore, the averaging problem in cosmology is of considerable importance for the cor- 
rect interpretation of cosmological data. The correct governing equations on cosmological 
scales are obtained by averaging the Einstein field equations (EFE) of GR (plus a theory 
of photon propagation; i.e., information on what trajectories actual particles follow). By 
assuming spatial homogeneity and isotropy on the largest scales, the inhomogeneities affect 
the dynamics though correction (backreaction) terms, which can lead to behaviour qualita- 
tively and quantitatively different from the FLRW models. It is necessary to use an exact 
covariant approach which gives a prescription for the correlation functions that emerge in an 
averaging of the full tensorial EFE. For example, in [IT] the macroscopic gravity equations 
were explicitly solved in a FLRW background geometry and it was found that the correlation 
tensor (backreaction) is of the form of a spatial curvature. 

Clearly, backreaction (averaging) effects are real, but their relative importance must 
be determined. In the FLRW plus perturbations approach, the (backreaction) effects are 
assumed small (and are assumed to stay small during the evolution of the universe). But it 
is also possible that averaging effects are not small (i.e., perturbation theory cannot be used 
to estimate the effects, and real inhomogeneous effects must be included). The Wilkinson 
Microwave Anisotropy Probe (WMAP) [2J, together with SNIa data in ACDM models [1], 
suggests a normalized spatial curvature Qk ~ 0.01—0.02 (i.e., of about a percent). Combining 
these observations with large scale structure observations then puts stringent limits on the 
curvature parameter in the context of adiabatic ACDM models; however, these data analyses 
are very model- and prior-dependent [8], and care is needed in the proper interpretation of 
the data. There is a heuristic argument that Qk ~ 10~ 3 — 10~ 2 [T2l [T3j . which is consistent 
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with CMB observations [13] and agrees with estimates for intrinsic curvature fluctuations 
using realistically modelled clusters and voids in a Swiss-cheese model [T5| [16] . It must be 
appreciated that such a value for O/., at the 1% level, is relatively large and may have a 
significant dynamical effect on the evolution of the universe; the correction terms change the 
interpretation of observations so that they need to be accounted for carefully to determine 
if a model may be consistent with cosmological data [HI [13] . 



2 Null geodesies 



All deductions about cosmology are based on light paths. Only the redshift and the energy 
flux of light arriving from a distant source are observed, rather than the expansion rate or 
the matter density. It is often assumed that intervening inhomogeneities average out. How- 
ever, inhomogeneities affect curved null geodesies [121 H5] and can drastically alter observed 
distances when they are a sizable fraction of the curvature radius. In the real universe, voids 
occupy a much larger region as compared to structures [17], hence light preferentially travels 
much more through underdense regions and the effects of inhomogeneities on luminosity 
distance are likely to be significant. 

The effect of averaging null geodesies in inhomogeneous models was discussed in [19]. 
GR is treated as a microscopic (classical) theory. Real photons travel on null geodesies 
in the microscopic geometry. However, because all observations are of finite resolution, 
observations necessarily involve averages of measured quantities. Therefore, in interpreting 
real observations, it is necessary to model properties of (not only a single photon but of) a 
'narrow' beam or bundle of photons (i.e., a local congruence of null geodesies). From the 
geometric optics approximation we can obtain the optical scalar (Dyer-Roeder) equations 
that govern the propagation of the local shearing and expansion (of the cross-sectional area 
of the beam) with respect to the affine parameter along the congruence due to Ricci focussing 
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and Weyl tidal focussing [TS]. Since the nonlinear optical scalar equations require integration 
along the beam, the optics for a lumpy distribution does not average and there may be 
important resulting effects. Therefore, it is important to study the effect of averaging on a 
beam of photons in the optical limit. 

The motion of photons in an averaged geometry (and the resulting effect on cosmological 
observations) was discussed in [19]. Assuming GR to be a microscopic theory on small scales, 
with local metric field g (the micro-geometry) and matter fields, a photon follows a null 
geodesic k in the local geometry. After averaging (using a covariant averaging scheme), we 
obtain a smoothed out macroscopic geometry (with macroscopic metric (g)) and macroscopic 
matter fields, valid on larger scales. But what trajectories do photons follow in the macro- 
geometry? In general, the "averaged" vector (k) need not be null, need not be geodesic 
(and even if it is, need not be affmely parametrized) in the macro-geometry. This clearly 
affects cosmological observations. In the case of radial geodesies in a sphericallly symmetric 
geometry [19] . it was found that there is an effect due to the non-affine parameterization of 
the null geodesies in the averaged geometry, together with additional small corrections due 
to the fact that the trajectories are not exactly null geodesies, the correlation tensor is not 
precisely due to a spatial curvature and since a beam of photons will experience expansion 
and shearing during its evolution [18J. Each of these effects are typically of order of about 
1%, which can add up and perhaps produce a significant observational effect. 

Similar issues have been discussed recently from a different point of view. It was con- 
jectured that in a statistically homogeneous and isotropic dust universe, light propagation 
can be treated in terms of the overall geometry (meaning the average expansion rate and 
average spatial curvature) if the structures are realistically small and the observer is not in 
a special location [12] . The relationship between the expansion rate, the redshift and the 
distance scale in such a universe containing non-linear structures was further investigated in 
[2D] , and it was shown that light propagation can be expressed in terms of averaged geomet- 
rical quantities, up to a term related to the null geodesic shear. In general, the null shear 
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is not negligible, and thus the Dyer-Roeder equations do not correctly describe the effect of 
clumping. Instead, the redshift and the distance are determined by the average expansion 
rate, the matter density today and the null geodesic shear. This implies that a clumpy model 
can be consistent with the observed position of the CMB acoustic peaks even when there is 
significant spatial curvature [2J, provided that the expansion history is sufficiently close to 
the spatially flat ACDM model [20]. 

In addition, the propagation of photons in a particular (toy) Swiss-cheese model, where 
the cheese consists of a spatially flat, matter only FLRW solution and the holes are con- 
structed out of a LTB solution of the EFE, and the phenomenological effects of large-scale 
nonlinear inhomogeneities on observables such as the luminosity-distance-redshift relation, 
were discussed in [9]. Following a fitting procedure based on light-cone averages, it was 
found that the expansion scalar is unaffected by the inhomogeneities (essentially due to the 
spherical symmetry of the model), but the light-cone average of the density as a function 
of redshift (which is not affected by the spherical symmetry) is affected by inhomogeneities. 
(The effect arises because, as the universe evolves, a photon spends more and more time 
in the (large) voids than in the (thin) high-density structures.) The phenomenological ho- 
mogeneous model describing the light-cone average of the density is similar to the ACDM 
concordance model, and behaves as if it has a dark-energy component. 

Clearly averaging can have an important effect on photon propagation and hence ob- 
servations. Ideally, we would like to describe the physical effects of averaging in terms of 
observational quantities, focussing on the behaviour of light. In particular, spatial averages 
at constant time are not directly related to observations; it is observationally more mean- 
ingful to consider light-cone averages of quantities. Consequently, we wish to attempt to 
describe the evolution equations on the null cone, focussing on observational quantities. Let 
us first consider the formulation of the EFE on the null cone. 
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3 EFEs in Observational Coordinates 



The cosmological data representing galaxy redshifts, observer area distances and galaxy 
number counts as functions of redshift are given, not on a space-like surface of constant time, 
but rather on our past light cone C~(po), which is centered at our observational position po 
"here and now" on our world line C. By using observational coordinates (OC) [21], the EFE 
can be formulated in a way which reflects both the geodesic flow of the cosmological fluid 
and the null geometry of C~(po), along which all of the observational information reaches 
us. In this formulation the EFE split naturally into a set of equations which can be solved 
on C~(po), that is on our past light cone, and a second set which evolves these solutions off 
C~(j?o) to other light cones into the past or into the future. The solution to the first set is 
directly determined from the data, and those solutions constitute the "initial conditions" for 
the solution of the second set. 

For example, the OC x % = {w, y, 8, 0} for a spherically symmetric metric are centered 
on the observer's world line C and defined in the following way: (i) w is constant on each 
past light cone along C; each observational time coordinate w = constant specifies a past 
light cone along C and our past light cone is designated as w = wq. (ii) y is the null 
radial coordinate, which measures distance down the null geodesies with afline parameter v, 
generating each past light cone centered on C, so that y increases as one moves down a past 
light cone away from C. (iii) 9 and (f) are usual the latitude and longitude of observation, 
respectively. In OC the spherically symmetric metric takes the form: 

ds 2 = -A(w, yfdw 2 + 2A{w, y)B(w, y)dwdy + C(w, y) 2 dtt 2 , (3.1) 

where the comoving fluid 4- velocity is u a = A -1 ^. There is remaining coordinate freedom 
which preserves the observational form of the metric. It is also important to specify the 
central conditions for the metric variables A(w, y), B(w, y) and C(w, y) in equation (13. ip (as 
they approach y = 0). The FLRW metric is obtained in OC by effectively specifying A = 
B = a(n),C = a(n)k~ 2 sin 2 (ky) in conformal time n = w — y. It was shown in [22] how to 
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construct flat dust-filled A ^ FLRW cosmological models from FLRW cosmological data 
on our past light cone, by integrating the exact spherically symmetric EFE in OC first on 
the light cone (integrating, e.g., ^-) and then off of it (integrating, e.g., ^). 

The basic observational quantities on C are then defined as follows: 

(i) The redshift z at time Wq on C for a comoving source a null radial distance y down 
C~(p ) is given by 1 + z = A(w Q ,0)/A(w ,y). 

(ii) The luminosity distance d^ [21], measured at time wq on C for a source at a null 
radial distance y, is given by di = (1 + z) 2 C{wq, y). 

(iv) The number of galaxies counted by a central observer out to a null radial distance y is 
given by N(y) = An f$ fi(wo,y)m~ 1 B(wo,y)C(wo,y) 2 dy, where /i is the mass-energy density 
and m is the average galaxy mass. 

Data is given on the null cone (in a suitably averaged form). In principle, we need 
to construct the metric on the null cone. For example, in the spherically symmetric case 
averaged quantities are given as simple integrations with respect to y (or redshift); matter 
terms, such as N(y), are already expressed as averaged quantities. To obtain the averaged 
geometry on the null cone (i.e., in terms of averaged metric functions (A), (B) and (C)), we 
still need to average the EFE on the null cone. Thus we still need to be able to average the 
tensorial EFE using some covariant averaging scheme, and then separate the equations into 
evolution equations on the null cone (in terms of (A), (B) and (C) and their derivatives), 
and constraints and evolution equations off the null cone. Therefore, we need to project the 
suitably covariantly averaged full EFE onto the null cone to compare with data averaged on 
the null cone. 

Thus, in this approach we are no further ahead, unless we can covariantly 'average' just 
an appropriate subset of the evolution equations on the null cone. Again we note that it is 
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not really the geometry that needs to be averaged, since it is not the geometry per se that 
is observed (or appears in the observational relations). That is, we observe the averaged 
matter on the null cone; if we average the geometry both on the null cone, and inside (and 
outside) of the null cone, we are averaging quantities that are not directly observable and 
may not play any role in determining the properties of the averaged matter on the null cone 
(either directly via the energy momentum tensor or the Ricci tensor through the EFE, or 
indirectly through the effects on the motion of matter in response to the geometry). Thus 
we would like to be able to focus on the effects of the averaged geometry that only affect 
cosmological observations. 

Suppose that it is possible to covariantly 'average' just a subset of the evolution equations 
on the null cone, focussing on relevant quantities (and not necessarily the complete averaged 
spacetime geometry). Perhaps we can focus just on scalar equations (and not the full set 
of EFE) in the spirit of the Buchert approach. In this approach to the averaging problem 
(for irrotational dust) only scalar quantities are averaged, yielding the averaged Hamilto- 
nian constraint (or generalized Friedmann equation), the averaged Raychaudhuri equation 
(plus an integrability condition) [15]. The Buchert approach is 'heuristic' since the Buchert 
equations are not closed. Indeed, since only scalar quantities are averaged not all of the 
EFE have been averaged in Buchert's approach, and consequently any solutions for these 
equations may not be consistent with the full set of EFE. 

4 Scalar Equations 

Let k a represent the tangent vector of the bundle of affmely parametrized null geodesies 
(optical rays) with scalar expansion (9 = \k a . a ) and shear (<r) [18]. The null version of the 
Raychaudhuri equation is then [23] 





D e = -d 2 -d 2 - 



—R a bk k , 



(4.2) 
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where DO = 9^ a k a is a directional derivative along the null geodesic (and can be expressed in 
terms of the affine parameter, v, or the redshift) and a 2 is a scalar. In the Newman-Penrose 
formalism, this is effectively eqn. (7.21a) in [23] and constitutes the first of the optical scalar 
equations (there is also an equation for D(a 2 )). 

Let us assume that the matter is dust: 

T ab = pu a u b . (4.3) 

We can decompose the fluid four- velocity as u a = k a + v a , where V = k a v a = k a u a (v a v a = 
— 1 — 2V) and is related to the redshift by 1 + z = V/Vq. We can then rewrite eqn. (14. 2p as: 

D6 = -6 2 -a 2 -\pV 2 . (4.4) 

The projected part of the conservation equation, T ab )b = 0, becomes: 

Dp = -p0-A, (4.5) 

where A = p^v b + pv ,b . We can write an "effective" first integral in the form 

l -P = pV 2 + B, (4.6) 

where D(B) is obtained from eqns. (j4.4p and (14.51) . 

Equation (14.41) concerns the propagation of 9 down the null cone. Let us average this 
scalar equation along the null cone: 

^(6} = -(e) 2 -yl(p}(l + z) 2 + Q, (4.7) 

where Q = (8) 2 — (6 2 ) — (B 2 ). We can write an averaged version of ( 14.61) in the form 

\{0f = V 2 (p)(l + zf -\r~\q, (4.8) 

where R is defined through this equation. 
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Equation (j4.7p is an important result; it is a scalar equation which formally represents 
how the averaged expansion varies along the null cone. It (together with an equation like 
(14. 8p ) represents a null version of the Buchert equations. It suffers the same disadvantages 
as the Buchert equations in that not all of the appropriate EFE have been averaged and the 
system is not closed. However, it does have the advantage of relating physically observed 
qualities (such as the density averaged on the null cone), and any assumptions to close the 
system now relate correlations on the null cone that may be physically better motivated. 

Alternatively, we can define a characteristic length scale I of the average area behaviour 
of the geodesies by: 

2 dv 
Eqn. (14. 4ft can then be written in the form 

d 2 t ^ 9 
= -20 2 - P V 2 . 

dv z 



4.0.1 FLRW models 



The spatially homogeneous and isotropic FLRW metric can be written in the form: 

d s 2 = -dt 2 + a 2 (t) [dr 2 /(l - kr 2 ) + r 2 d6 2 + r 2 sin 2 6d(j) 2 ] , (4.9) 

where a(t) is the cosmic scale factor and k is the curvature of 3-dimensional space (e.g., k = 
corresponds to a spatially flat universe). The wavelength A of photons moving through the 
universe scale with a(t), and the redshift of light emitted from a distant source at time t em 
is defined by 1 + z = A b s /A em = l/a(t em ); thus dt = —dz/H(z)(l + z), where H = a/a is 
the Hubble parameter and an overdot denotes a time derivative. For an object of intrinsic 
luminosity L, the measured energy flux F defines the luminosity distance dh to the object, 



where di(z) = J L/AnF = (1 + z)r(z) , and r{z) is the comoving distance to an object at 
redshift z. The key equations of cosmology are the Friedmann and Raychaudhuri equations: 

cy (a\ 2 1 k A n . 
\aj 6 a z 6 

11 



- = -g </>+*) +3 (*■»> 

where p is the total energy density of the universe, p is the total pressure, and A is the 
cosmo logical constant. For each matter component i, the separate conservation of energy is 
expressed by d(a 3 pi) = —pida 3 . The deceleration parameter, q(z), is defined as q(z) = —^jj2- 

As an illustration, let us consider equation (14.41) in the case of a pressure-free, flat (k = 0) 
FLRW model with metric (14. 9p . where 1 + z — I /a. For 

k a = Q,1,0,0), (4.12) 

we have by direct calculation that 

9 = h a , b g ab = ~{a--), (4.13) 
2 a z r 

and 

a 2 = \k {a , b) k^ - d 2 = 0, (4.14) 

whence equation (14.31) yields 

ad — a 2 = —-pa 2 , (4-15) 
as expected (and corresponding to eqns. (14.101 - I4TTTT) ). 
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